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I. INTRODUCTION 



A wormhole is a solution for the Einstein's gravitational field equations that represents a hypothetical tunnel 
through spacetime. The detailed calculations about them are not simple, nevertheless with the techniques taught in 
a standard calculus course pQ, it is not difficult to produce the embedded diagram of a wormhole based on geometry 
principles [2]. 

Since the discovery of the Einstein equations, several features were addressed about these solutions, such as Einstein- 
Rosen bridges [3]; Wheeler [4] analyzed these kind of solutions and called them wormholes. The most important and 
popular contribution to this solution was given by Morris and Thorne [5] , they were the first who made a complete and 
detailed investigation in this field. Wormholes are solutions of the Einstein field equations which represent two connect 
universes or a connection between two distant regions of universe. The most important solutions, which represent 
a static spherically symmetric spacetime possessing two asymptotically flat regions are: the Lorentzian (pseudo- 
Riemannian manifold) and the Euclidian (real Riemannian manifold) wormholes [6 . Both of them are interesting, 
but the real physics take place in the Lorentzian type. For realistic models of traversable wormhole, one should address 
additional engineering issues such as tidal effects. People are sceptic about them because they contradict some energy 
conditions, for example, they need matter with negative energy density as a source [5]. Visser claim in his book [BJ: 
A wormhole is any compact region of spacetime with a topologically simple boundary but a topologically nontrivial 
interior. The existence of them imply violations to the null energy condition as well as causality among others strange 
features. In quantum mechanics, there are many systems that violate this condition, but understanding the relation 
between gravity theories and quantum mechanics is in debate in these days. However, in this paper we are assuming 
that wormholes exist and we analyze the stability properties of the most simple of them. This is important because 
recently it has been increased the investigation about its behavior. A full revision about these solutions after the 
Visser 's book [B] was done by Lobo [7]. 

About lineal perturbation theory in wormhole solutions, some of the most important works have been done by Frolov 
and Novikov [5], they analyzed them when interacting with an external electromagnetic or gravitational field. Kar 
et al [9 ], studied the propagation in the wormhole geometry of different types of perturbations but not gravitational 
perturbations. They studied the reflection and transmission of massless scalar waves in presence of an ultrastatic 
wormhole. Perez and Hibberd [ID] , studied the same effect about the electromagnetic waves using numerical tools. 
A more detailed work about perturbation theory for a ultrastatic wormhole using Teukolsky formalism, was done 
for Torres del Castillo and Lopez-Ortega They analyzed the scalar, neutrino and electromagnetic perturbations 
using adjoint operators method for a ultrastatic wormhole spacetime. With this method Torres del Castillo et al found 
similar results in the electromagnetic case than Perez and Hibberd, which give us a very good signal of the confidence 
of the Debye method. They obtained the transmission coefficient of a Schrodinger type equation, establishing the 
relationship between this coefficient and the amount of radiation transmitted through the wormhole. Finally, Kim 
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[12) studied the perturbations of the scalar, electromagnetic and gravitational waves under the traversable Lorentzian 
wormhole geometry; Kim found a Regge- Wheeler type equation for gravitational perturbations. Using non-linear 
perturbations, the instability of wormholes has been shown in [13j but only numerically. 

In this paper we use Ncwman-Penrose formalism to analyze the gravitational wave perturbations in the geometry of 
an ultrastatic worhmholc using Debye potentials. Also, we studied the propagation of gravitational fields showing that 
the radial part of the perturbations can be expressed in terms of Schrodinger type equation. In sec II the wormhole 
spacetime is described. In sec. Ill the ultrastatic wormhole geometry is described in terms of a null tetrad. In sec. IV 
we find a Schrodinger type equation for the perturbed Einstein equations. We analyze the reflection and transmission 
coefficients of the gravitational waves as well as the conservation of the energy. Finally conclusions are drawn in sec. 
V. 



II. WORMHOLE SPACE-TIME 



One solution to the Einstein equations which represent static spherically symmetric spacetime possessing two 
asymptotically flat regions is a Lorentzian wormhole, its metric has the form 

rlr 2 

ds 2 = e 2 ^dt 2 - - v 2 (d0 2 + sin 2 ddep 2 ) , (1) 

where <I>(r) usually is called the redshift function, because it is related to the gravitational redshift. The b(r) function is 
related to the shape of the wormhole using embedding diagrams. They satisfy the conditions: i) ^p- < 1, b(r = ro) = 

bo', ii)r — > oo, — > 0; iii)<fr{r) always finite. The radial coordinate has a range that increases from a minimum 
value at ro (where b(ro) = ro), corresponding to the wormhole throat, to infinity. As r — *■ oo, b(r) approaches 2M 
which is defined as the mass of the wormhole. 

The Einstein equations imply that such a metric exists if and only if the energy density and the radial tension in 
the frame of a static observer satisfies the relation p — f < (p denotes the energy density and f denotes the radial 
tension), at least for a small region that includes the throat of the wormhole f denotes the negative of the radial 
pressure. 

For a traversable wormhole it must have no horizons which implies that must be finite everywhere. To keep 
the radial tides being zero it is sufficient to set <&(r) equal to zero, resulting in an ultrastatic wormhole. In this case, 

we take a particular shape function b(r) that satisfies the wormhole requirements b(r) = -S, where bo is a constant. 
Therefore, the metric that we will analyze is the following 

ds 2 = dt 2 - - v 2 (d9 2 + sin 2 Stiff) . (2) 

The proper radial distance is defined by, 

r dr 

«r) = ± -j^=, (3) 
Jro . i _ Mr) 



measured by a static observer; this distance must be well behaved in any region of the spacetime. As is well known, 
the metric Q can be written in terms of the proper radial distance by 

ds 2 = dt 2 - dl 2 - r(l) 2 (d9 2 + sin 2 6d<p 2 ) , (4) 

where r 2 (l) — l 2 + bg, in the general case this metric represents a spherically symmetric static wormhole provided that 
(i) I ~ + oo , r 2 (l) oc I 2 ; (ii) I — > 0, r(l) — * bo; (in) <&(/) always finite. These conditions arc fulfilled if we take expressions 
for the ultrastatic wormhole. The proper distance decreases from I = +oo in the upper universe, to I — at the throat, 
and then from zero to / = — oo in the lower universe. We have changed the signature of the metric according to |14j 
to use the Newman-Penrose formalism. This method let us to decouple the perturbed Einstein equations determined 
by one angular equation and another radial Schrodinger- type differential equation. The potential associated to this 
latter equation gives us information about the stability of the system and their solution give us the energy fluxes. In 
this paper we study the gravitational perturbations of the wormhole metric following the method used in |llj . This 
method consists in the use of vector potential perturbations of a solution of the Einstein equations expressed in terms 
of a scalar (Debye) potential. 
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In the ultrastatic wormhole metric, the gravitational perturbations are analyzed as a superposition of different 
modes with a i-time and azimuthal </>angular dependence specified by 



e -i{u>t-m<p) ) (5) 

where m is an integer and w is a real positive constant. 

For a description of the wormhole space-time solution Q in a Newman-Penrose formalism |14j . the metric is 
described by 

9nv = 2[l<jJh,) - mQjjn v )], (6) 
the bar denotes complex conjugated. The null tetrad l,n,m and m is defined by the following identities, 

l-n=l and m-m=— 1. (7) 

So, we first need to define a null tetrad |15j according to the wormhole metric, which is satisfy by 

P=i(l, 1,0,0), 

n"=i(l, -1,0,0), 

= ~= — - — — (0, 0,1, i esc 9) , 

= — — — —(0,0,1, -i esc 6»), (8) 

V2(P + fog) 1 / 2 ' 

or equivalently we can define the directional derivatives, 



D = 


= 




f di), 


A = 








5 = 






1 




+62)1/2 


5 = 


m M <3 M = 




1 




+ 62)V2 



(dg + icsc9d v ) , 

(9e — i esc 6 d^p) . (9) 

The I?, A, <5 and 8, when acting on functions with a £ dependence of the form e~ luJt , are still just derivatives of I. 
The only nonvanishing component of the Weyl scalar in the background spacetime is 

*2 = (10) 

3(/2+62)2 

which means that the wormhole metric solution is type D in Petrov classification [14 . In order to decouple the 
perturbed Einstein equations (radial and angular part), the Newman-Penrose formalism is applicable only in type D 
metrics. Type D space have two double principal null directions which define "radially" ingoing and outgoing null 
congruences near the object which is the source of the field. 

In the wormhole metric, the spin coefficients associated to this geometry are given by 

h> = (7=K = \ = e = "f = T = TT = 0, 

I 

p = fj, 



V2(P + 6§) • 



P=-a=^ 2¥ ^ = S^0, (11) 
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where the only nonvanishing Ricci scalars are 



$oo = - n „, ,°„„ , $11 = 



2(Z 2 + 6 2 ) 2 ' " 2(/ 2 + 6g) 2 ' 

ft 2 ft 2 
22 2(Z 2 + 6 2 ) 2 ' 12(Z 2 + fog) 2 ■ [ 1 

Finally we can see the next relations between the nonvanishing component of the Weyl spinor and the Ricci scalars, 

2 4 2 

*2 = g*oo = -3*11 = 3*22 = 4A, 

here A is proportional to the curvature scalar R — g^ v R^ v ■ 

Now we already have all the ingredients for use it in the adjoint operators method to find the gravitational pertur- 
bations of the ultrastatic wormhole metric Q . 

III. GRAVITATIONAL PERTURBATIONS 

In this section, we will use the perturbed Einstein equations to obtain the Teukolsky equation that allow us to 
analyze the stability of the wormhole solution through the radial equation and its potential. 

In the Newman-Penrose formalism there are six equations - four Bianchi identities and two Ricci identities - which 
are linear and homogeneous in the quantities which vanish identically in the background, (see §29 of |14j , (321(a)), 
(321(e)), (321(d)), (321(h))): 

(6 - 4a + tt)* - (D - 2e - 4p)*i - (3* 2 - 2$ u ) re = (8 + 7f - 2a - 2/3)<J> 00 - (D - 2e - 2p)$ i - k$ 2 + 2<7$i , 

(13) 

(A - 4 7 + p)* - (8 - 2/3 - 4r)*i - (3* 2 + 2$ u ) a = (8 + 2fr - 2/3)<J> i - (D - p - 2e + 2e)$ 02 - 2k<J> 12 - A$ 00 , 

(14) 

(D + 4e - p)* 4 - (5 + 4tt + 2a)* 3 + (3* 2 + 2$n) A = (5 + 2a- 2r)$ 2 i - (A + p, + 2 7 - 2 7 )$ 02 + 2^<J>i + d$l%) 
(8 + 4/3 - r)* 4 - (A + 2 7 + 4p)* 3 + (3* 2 - 2$ u ) z/ = (5 - t + 2a + 2/3)$ 22 - (A + 2p + 2 7 )$ 2i - z?<J> 20 - 2A<£1§) 

and, on the other hand, the equations (310(b)) and (310(j)) from [H] are respectively: 

*o + (5 - 3/3 - a - t + W)k - (D - 3e + e - p - p)a = 0, (17) 
* 4 + (A + p + /I + 3 7 - 7) A - (8 + 3a + /3 + tt - f )i/ = 0. (18) 

These equations are linearized in the sense that the Weyl scalars ^4 an d the spin coefficients k, a, A, 

and v as perturbations, are taken up to first order only. We shall arrive to equations that involve ^0 an d *4 since 
they are gauge invariant in a linear perturbation theory. This components are the most significant in the gravitational 
radiation [16] . 

Since the ultrastatic wormhole metric is type D, then the non-perturbed quantities K i cr, A and p vanish, so 



we have from equations ( 13 1-(|16[), 

(5 - 4a)*f - (D - 4p)*f - (3* 2 - 2$ u ) n B = - (D - 2p) + {8 -2a- 2/3)$f + {8 -2a- 2/3) B $ 00 , (19) 

(A + p)^ B -(8- 2/3)*f - (3a B * 2 + 2$ u ) <r B = (8 - 2f3)$g - (D - p)** - A s $ o, (20) 

(8 + 2a)*f -(D- p)¥f - (3* 2 + 2$u) A s = (A + p) $f - (5 + 2a)*?! - <r S $ 22 , (21) 

(A + 4/x)*f - (5 + 40)¥f - (3* 2 - 2*n) v B = (A + 2/2)$^ - (J + 2a + 2/3)$f 2 - (* + 2a + 2 / 3) B $ 22 , (22) 



here we are used (111 and the superscript B denotes the perturbed quantities. 



The perturbation of the Bianchi equations (17) and (18 1 are, respectively 



{D-4p-p)^ 2 a B -{8-Zp-a)y 2 n B -^ B ^ 2 = [2p$n - (A + p) $ o - 2DA] a B (23) 
(J+3a + ^)* 2 ^ s -(A + 4p + 17)* 2 A B -*f* 2 = [2 M $ 11 + (D-p)$ 22 + 2AA]A B . (24) 
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where we are used the equations (321(b)) and (321(g)) from [T3] . 



If we operate (8 — 3/3 — a) on equation ( |19[ ) and with (D — Ap — p) on equation ( 20 1 , and subtract one equation from 
other. Finally we use the relation and the identity gB] [(£> -4p-p) (8 -2/9) - (6 - 3/3 + a) (D - Ap)} = to 
vanished the terms in , we obtain the Teukolsky equation for 



[(D - Ap - p) (A - 4 7 + p) - {5 - 3/3 - 5) (J - 4a) - 3$ 2 ] $f = T + T 0a , 



(25) 



where 



To = (8 - 3/3 - a) [(£> - 2p) - (5 - 2a - 2/3)** ] + (£> - 4p - p) [(5 - 2/3)$* - (U - p)<] , (26) 

(T is equivalent to (2.13) in (55]) and the extra term 

T 0a = [2 (D - p - p) $ n - 3 (A + p) $ o - GDA] a B + [2(5 - 3/3 - a)$ u ] k b - [(£> -Ap-p) $ o] A s 
-(<5 - 3/3 - a){6 -2a- 2/3) B $ o, 

represent the source when the perturbed differential operator 8 and the perturbed spin coefficients <r, ft, A, a,/3 are 
non- vanished [19 . 

Similarly to *S?q , we obtain the Teukolsky equation for "J J 3 using the relation (24 1 and the relation 
[(A + Ap + ft) (6 + 2a) - (8 + (3 + 3a) (A + Ap)] = to eliminate *f , we obtain 

[(A + Ap + ft) (D - p) - (5 + + 3a) (8 + 4/3) + 3$ 2 ] f = T 4 + T 4a , (27) 

where 

T 4 = (5 + /3 + 3a) [(A - 2p) - (5 + 2a + 2/5)$f 2 ] + (A + 4/z + Jt) [(5 - 2a)0>f 1 - (A - ^)$f ] , (28) 
(T 4 is equivalent to (2.15) in [55]) and the extra term 

T 4a = [3 (A - p) $ 22 - 2 (A + p + p) $ u + 6AA] A B -2(a+/?+3a)*iii/ B +(A + 4 M + /2) <J>2 2 cf S -(^+3a)(<$+2a+2/3) s $ 2 2, 

represent the source when the perturbed differential operator 8 and the perturbed spin coefficients \,p,a,a,/3 are 
non- vanished. 



The Ricci tensor terms on the right-hand sides of equations (26 1 and (28 1 given by the Einstein field equations are, 



$00 = --RvvW = ^nT^lu = AirTu, 



(29) 



and so on, where R^ v is the Ricci tensor and T„ w the stress-energy tensor. 

Notice that the non-zero projections are functions of r and t only. Thus substituting the projections into the RHS 
of (25) and (27 1, the only Ricci terms which remain are $00 j $11, an d < &22, that are related by Tu, T nn and T mm . 



However these operators are now purely angular operators, and since the projections are functions of (t,r), the result 
is null (zero). The expressions for T and T 4 vanished. By the same way, from the relations given in the appendix 
|19j . the expressions for Tq u and T 4a are zero. In this way, we conclude that radially infalling dust into the wormholc 
does not produce perturbations on the perturbed Weyl scalars and $f . 

The Newman-Penrose formalism, allow us to reduce the Einstein equations just to four perturbed equations whose 
solutions are given by the Weyl components: 



V2 



V2., lB 



(30) 



From each of these decoupled equations (25l-(27l and applying the adjoint operators method [20], one obtains a 
local expression for the solution of the Einstein equations in terms of a scalar potential. The equations are 



[(A + ft) (D + 3 P ) - (8 + 3a + /3) (8 + 4/3) - 3$ 2 ] ^fX = °> 



[(£> - p) (A - 3p) - {5 - a - 3/3) (5 - a + /3 + 2/3) - 3* 2 ] if> 



(G) 

(4) _ 
(G) - 



0, 



(31) 
(32) 



where tp(Q\ and iP)q\ are the scalar potentials associated with the decoupled equations for *S>q and Hf^ in the form 



',(4) 



y (G) 



^ 0) - $f and ^ 4) 



scalar potential tpo 20J , which satisfies the equation 



p 4 $ 4 3 . We can write the equations (31 (-( 32 ) as one single equation in terms of the complex 



[(A + J1)(D + 3p) -(S + 3a + /3)(<5 + 4/3) - 3$ 2 ]^ G = 0, 



(33) 
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substituting ipo = ipQ and ipa = m ^ ne equation (33), we are able to obtain the equations (31 1 and (32 1 



respectively. 

As it shown in [HJ-J53], the outgoing and ingoing metric perturbations of an algebraically special solution of the 
Einstein vacuum field equations in a frame such that = ^1 = (or equivalently k = a = 0), are given by 



: 2{l ll l„[(6 + 3/3 + a - t)(S + 4/3 + 3r) - \{D + 4e + 3p)] + m^m^D + 3e - e - p)(D + 4e + Bp) 
l^m v) [{D + 3e + e - p + p)(<J + 4/3 + 3t) + (<$ + 3/3 - 5 - r - tt)(£> + 4e + 3/j)]}^ ( ( §j + c.c, 



(34) 



2{— n„n M (<5 — a + [3 — n)(5 — a + 3/3 — 67r) — m M m iy (A — 37 + 7 — /-t)(A — 47 — 3/i) 



+ [(A - 37 - 7 + /Z - /z)((5 - a + 3/3 - 6tt - r) + (<5 - a + 3/3 - n + t)(A + 47 + 3p)]}^ + c.c, (35) 



respectively. 



In order to obtain a radial equation, we look for separable solutions of equations (31 I and (32 1 in the form 



^,=e- fai _ 2 r;. m (^)l( )(!), 
^ ( g ) = e- <w *2^ m (e >¥ »)xW(0 l 



(36) 
(37) 



where the dependence in tp is implicit in the spin weighted spherical harmonics 2Y jm{6, f) =-2 Yj m ( 9, ip). 

If we substitute the directional derivatives djjj, the Weyl scalar (10 1 and the spin coefficients pTj) into equations 



(31 1 and (32), we find that the functions X^(T) and X^(l) must satisfy the equations, 
21 



A, - 



2iujl 



5« 



(i + 2)(j-l) 



e-- t _ 2 F J - m (e,^)X(°)(0=0, 



(38) 



2/ 



A_ 



2iujl 



5bl 



where we have used the definitions 

d 



u+w i y\ e -i»t 2Yjm{dj(p)x w {l)=0i 



dl 



A = A + A_ = -jp+u 



In terms of a new function i?(/)( 2+s ) (where s is —2 and 2 for and X( 4 )(Z) respectively) defined by 

a(o (2+s) = z(o (2+s) a 2 + &o) 1/2 : 



we obtain a master Schrodinger-type equation from (|38| and (39 1 

2iwls (2Z 2 + 4fc§) (j + 2)(j-l) 



P + bl {P + biy 



p + bl 



z^ +s \i) = 0. 



(39) 



(40) 



(41) 



(42) 



Thus, and A( 4 )(Z) satisfy complex-conjugate equations (provided that u> is real, which is the only case to 

be considered here). 

Using spin weighted spherical harmonics, the dependence in jm is represented by the term ^'"ffi^T 1 ^ in the last 



equation, and since both equations ( 38 1-( 39 1 satisfy conjugated complex 



J(p=Q 



(43) 



then, equation (42) depends only on the radial coordinate I. 



Then, equation (42 1 can be transformed into a one-dimensional wave equation of the form, 

A 2 Z {2+s \l) = V{l)Z^ s \l). 



(44) 
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As we can see, V(l) is a function with two kind of potentials, one real and one imaginary 

2Juds (2/ 2 + 4bg) (3+2KJ-1) 

lj ~ ; 2 + 6g (J 2 + &§) 2 i 2 + &§ ' 1 j 

the real potential (for j > 0) is of short range, have a very well behavior and vanish at infinity. 

Then, the perturbations to the wormhole geometry is governed by a single Schrodinger-type one-dimensional wave 
equation 

d ^ + u 2 Z = V J Z, (46) 



where the potential Vj depends on the angular momentum j of the perturbations. Hence, equation (46 1 admits 
spatially well behaved solutions when u> is positive, being ui 2 always positive. Solutions of this equation with negative 
ui 2 correspond to exponentially growing and decaying modes which imply instability. 



We can rewrite equation ( 46 ) as 



oy , Z=-^- + V j Z, (47) 



if the operator on the right-hand side of (47 1 is positive, which is the case if Vj is positive definite, it follows that u> 
is real and positive and there are no unstable modes (if u> = there might be modes with linear growth in t) . On the 
other hand, if Vj is negative the operator is not necessarily positive. If one can find an eigenfunction with negative 
eigenvalue it is complex and one has a solution that grows exponentially in time. In our case, the real part of Vj is pos- 
itive for j > 0, and hence does not admit any io 2 < 0, which tell us that the wormhole is stable in accordance with [2~4"] . 



Multiplying equation (47 1 by Z (complex conjugate of Z) and take the integral from — oo to oo respect to coordinate 
I, we have: 

/+oo r+oo /"-t-oo 

(ReV-oj 2 ) \Z\ 2 dl+ / {imV -u 2 )\Z\ 2 dl + |Z'| 2 d/ = 0, 

-oo J— oo J —oo 

because the potential have real and imaginary part. In the first case (first integral), we can observe that ui 2 have to 
be real. If the imaginary part is odd function respect to w, we can argue that Im (iuj) — Re {ui). Therefore, it excludes 
an exponential growth of any bounded solution of equation (42). From Fig. [I] we can deduce that the integral of the 



imaginary part of the potential is zero, since is an odd function independently whether s — 2 or s — —2, reason why 
it does not contribute to the energy integral. So, the imaginary part does not give us information about the stability. 



Let Z{1) denote the solutions of one-dimensional wave equation with a bounded short-range potential V{1) (45 1, 
with the following asymptotic behavior when I — > ±oo 

Z — ► F s e ±iul , 

hence, we can assume that when I — > +oo the functions Z^ and Z^ have the form 

~— iujl 

Z (0) _^ A ( ° )e — + B^Pe*" 1 , (48) 
Z W —> A^l 2 e~ iul + flWy, (49) 

these boundary conditions correspond to an incident wave of A^ amplitude from +oo giving rise to a reflected wave 
of amplitude in the Z^ wave; A^ and B^ for the Z^\ respectively. As we showed ZW and Z^ satisfy 



complex-conjugate equations (44 1 and their physical interpretation is the following: the Z^ wave correspond to the 
ingoing wave, i.e. the incident wave is approaching to the wormhole and the reflected wave is moving away from the 
wormhole; the Z^ correspond to the outgoing waves and in this case the wave moves in opposite directions to the 
wave. This behavior is analyzed only in the far zone. 
In the same way that in the potential barrier in quantum mechanics, when we have two (independent) solutions, 
in this case Z^ and Z (4 \ the Wronskian [Z°,Z 4 ] = Z 4 A + Z° - Z°A + Z 4 must be constant [Hj. So, we obtain a 



Potential 
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Real potential j-0 
Real potential j=1 
Real potential j=2 
Real potential j=4 
maginary potential 



FIG. 1: Real and imaginary potential of the Schrodinger equation associated to gravitational perturbations of the ultrastatic 
wormhole metric. 



relationship between the coefficients in the equations ( 48 1 y ( 49 ) . This latter will give us the information about the 
energy conservation between the incident and reflected wave. Evaluating the Wronskian when I — ► +oo, we obtain 

Q = cte., (50) 

-oo we impose the 



A 2 ^2 



R (0) R (4) 



where Q is a constant different of zero. Similarly on the other side of the throat when / 
boundary conditions, 



Z (0) 



D 



(0)< 



-Ztol 



I 2 



(51) 
(52) 



where correspond to a transmitted wave through the wormhole in the ingoing waves and D% in the outgoing 
wave, and computing the Wronskian of Z^ ' and Z( 4 ) in this limit we have 



(53) 



Because we have the same constant in both equations (50 1- (53) we can write the next relation 



A i0) A (4) 

A 2 A 2 



(54) 



The reflection and transmission coefficients for both potentials V(l) must satisfy 



ii?r 



A 2 A 2 



\T[ 



4(0) .(4) 
A 2 A 2 



(55) 



then, from equation (54) we have the standard relation for the energy conservation in the potential barrier in quantum 
mechanics 



|i?| 2 +|7f = l. 



(56) 
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IV. REFLECTION, TRANSMISSION AND ABSORPTION OF GRAVITATIONAL WAVES 



In section III, equations (31 l-(32) was used to obtain a decoupled equation (radial and angular) from the Teukolsky 
equation. From equation (54 1 we found the reflection and transmission coefficients (55). 

In this section we will to calculate the reflection and transmission coefficients in terms of the explicit form of the 
perturbations of the Weyl scalar distinguished with a superscript (B). To do that we calculate making use of (34 1- (35 1 
and the relation 



*ACDE = ^(A^chDE)R'S' + jh^AODE, 



1. 



?R'v7S' 



(57) 



f g\ ( B) 

the components "Jq and ty\ of the gravitational perturbations, where in this case the trace h is zero [25j . From 
this latter, the components obtained for the outgoing and ingoing waves are |23j : 



(B) _ -2Yjm{6, <j>) e" 



(B) 



j(j 2 -l)(j + 2) 

i(l 2 +b 2 ) 2 



d 4 t d 3 „ 2 dl 2 



^—iujt 

e 2 i j 



Y ]m (6,<f>)X (0 \l) 



dl 
Ibl 



H 



e iut X w (l) 2 Y j7 



(58) 
(59) 



(B) 



K^-1Kj + 2) 

4(/ 2 + 6g)2 

■iuji 



^ — ILOt 

e - 2 i j 



Y jm {e,4>)x^\i) 



Ibl 



662 



^(B) _ 2Y jm (8,<p)c 

T 4 A 



d*_ 
dl 4 



d?_ 

dl 3 



6w z 



,dP_ 
dl 2 



{ll + blY 3(/2 + 62)4 

d 



e iut X (i \l)_ 2 Y jm (e,, 



dl 



(60) 
(61) 



respectively. 

By means of the Teukolsky-Starobinsky identities, which in the present case take the form 



(D-3e + e- 7 p) (D + 2e-2e- bp) (D + 3e - e - 3p) (D + 4e - p) * 4 = 

(5 - 3/3 - a + 7n) (5 - 2a - 2/3 + 5vr) (5- 3a-/3 + 3tt) (5-Aa + tt) W - 3# 2 (/oA + pD - t5 - ttS)^ , 
(V + 7p + 37 - 7) ( V - 2 7 + 27 + bp) (V + 3/x - 37 + 7) (V - 4 7 + p) * = 

(5 + - 7t + 3a) (6 + 2(3+ 2a- 5r) ((5 + 3/3 - 3r + a) ((5 + 4/3 - r) * 4 - 3* 2 (pA + pD —tS — 7t5)W 4 , (62) 



we can find a relationship between the coefficients appearing in (|48|-(|49l 

A (f>) Ai) 

/In , /li 



(0) 



Considering the asymptotic expressions when / 

Z^(l) - ' 



2 ' 

DC 



4(0)' 
Ik 



\k=2 



/fc-4 



B. 



-iujl 



— iujl 



B. 



(4) 



E Sfe 



(0)' 



\k=2 



lk-4 



(4)' 



vfc=2 



\fc-2 



and when I 



/fc 



-iujl 



\k=2 



(4)~ 



\fc=2 



fc-4 



(63) 

(64) 
(65) 



and substituting Z(Z) in X(Z) (41 1 respectively, we can obtain de Weyl scalars ( 60 1-( 61 ) . 

Using the Teukolsky-Starobinsky identities (62) and the last asymptotic expressions, we obtain 



jU - 




-2)(j- 


d4 4) 


= 4^4°), 


j (3 - 




r2){j- 


1)^ 4) 


= 2(Bi 0) - 


jU - 




-2)(j- 


1)^ 4) 





?,2 R (0)u2 



(66) 
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As you can see, there are mixing between the coefficients A { 2 ] , A { 2 ] , B { 2 ] , B { 2 ] and B { 4 4) . So, now we will manipulate 
the coefficients in order to find the reflection and transmission coefficients which appear in the energy conservation 
equation (56). 



So, it is possible to find the coefficients A^?\ 



B 



(0) 



respectively, making use of (42 1. Thus, we find 



|(0) 



and Bi 1 for k > 3, in terms of A 



(°) 



A?> and B 2 



(4) 



A 



(o) 



(o) 



B 



\j(j + l)-6]A 



A 



(0) 
2 



2icu 

1) - 12]A 



(o) 

3 



(0) 



(o) _ 



[2 + !)]£. 
2mj 

-iCj + i)^ 0) - 

4iw 



Aiwb 2 B 2 0) 



(67) 



.4 



(4) ^ [ 7 -Q- + l)-2]4 4) 
2iw 

_ j(j+l)4 4) +4za;bg4 4) 



4 



M 4) = 



Bi 4) = 



(4) 



4zw 
2iu) 

[12-j(j + l]\B^-^Bl 



(68) 



Finally, with the equations ( 66 1-( 68 ) we obtain 

I (4) _ 



4wM 



4j(0) 



M 4) = 



M 4) = 



[i(i + i)(i + 2)(i-i)]- 



4w 2 



4w 4 L> 



4n(0) 



b'(j + l)0" + 2)(j - 1)]' 



(69) 



Substituting ( 69 1 into ( 54 ) it follows that 



|A 



l-D. 



(o), 2 ^ 6giO' + l)(j-l)0' + 2), B (o),a 



16w 6 



(70) 



from this last equation, we can find the expression for the incident, reflected and transmitted wave. 

The expression ( 70 ) was obtained from the formalism of the quantum mechanics although often it does not have 



tangible physical meaning. In the section "flux energy" this conservation will calculate directly from the point of 
view of Teukolsky equation with the wave functions that represent the gravitational waves, which is physically more 
significant. 



perturbations corresponding to the potential V(l) when I 



Expressing the radial functions X(l) of equations ( 58 1-( 59 1 in terms of ( 63 1-( 65 1 , the components of the gravitational 

oo are 





(B 



-2IWM), 



) j(j2-i)(j + 2) (o) g-M«-0 



(71) 
(72) 
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when I 



-oo, we have 



(B) 



4n(o)e -MW) 



-2Yj m (9, ' 



I 5 



(73) 
(74) 



as you can see the last component indicate us that there are no waves that leave the throat of the wormhole. 

A. ASYMPTOTIC BEHAVIOR 

In the outgoing gauge, since the potentials V(l) are of short range, Y(l) have the asymptotic behavior e ±lul both 
for I — > ±oo. If the function Y(l) have the asymptotic behavior e~ lbjl for I — ► oo, from equations ( 58 )-( 59 1 we have 



^ = O 



1 



l(l+a) / ' 



which means that for ingoing waves where a — 0, represents the ingoing fields. Similarly if Y{1) have the 
asymptotic behavior e lujl for I — > oo, 



= Q 



l(5-a) 



hence, at spatial infinity for outgoing waves where a = 4, tyf represent the outgoing fields. 



B. FLUX ENERGY 



As we know, the energy tensor of the gravitational field is not well defined, and in order to find the energy 
conservation, we use the formalism used by Chandrasekhar [14] ■ We combine the eq. (493)-(494) of the chapter 98(b) 
of |14j and eq. ( 71 1 to obtain the flux of energy per unit time due to the gravitational radiation coming from infinity, 



for the ingoing waves: 



dt 



^ lim 



(B) 



/dE?r\_u 

\ dt 



-in 



A. 



(0) 



(75) 



where the symbol () represent time-averaged. Similarly, the time-averaged flux per unit solid angle due to the outgoing 
gravitational wave (72 1, is 



d 2 E 9rav 
dtdil 

arj out 

dt 



1 



„ lim 

4:7TtO T— >oo 

bO' 2 - i)Q- + 2)| 

647TW 2 



D. 



(0) 



(76) 



Finally, using (73 1 we obtain that the flux of energy across the throat, is given by 

1 



dE grav 



dt 



( dK 
\ dt 



rav 
trans 



, lim 

47TW^ l—f—oo 

D. 



4 



4tt 



' ; (o) 2 

2 



Combining the equations (75 1- (77 1 through the conservation of energy flux 

arj trav° \ / aCj - 



arj out 

dt 



dt 



dt 



(77) 



(78) 



12 



we obtain the following expression: 

J(J + 1)(J-1)(J + 2)|^ 0) | 2 |Df| 2 

16^ 6 |4°)|2 |4°)|2 



= 1. (79) 



The energy conservation is fulfilled when the equations (79 1 and (70 1 are equal, i.e. 6q = 1 



V. CONCLUDING REMARKS 

Using the scalar Debye potentials for massless fields, the behavior of the gravitational perturbations of an ultrastatic 
wormhole are obtained in a similar way than for perturbations of black holes [2], [35] without sources. This work comes 
to complete the perturbation analysis made previously: electromagnetic, neutrinos and scalar fields perturbations [10 - 

We analyzed the potential associated to the Schrodinger-type equation that is obtained from the application of the 
Teukolsky formalism, and we concluded that the ultrastatic wormhole is stable under linear gravitational perturbations 
in agreement with |24j . 

In ref. [13 nonlinear perturbations was made and they found that the wormhole must be unstable under second 
order perturbations but they do not obtain any physical amount like energy conservation. We calculated the reflection 
and transmission coefficients from quantum mechanical point of view to verify the energy conservation. On the other 
hand, we calculated the energy conservation relation through energy fluxes using the Weyl components that represent 
the ingoing and outgoing gravitational waves. Comparing both energy conservation relations, we can conclude that 
energy conservation is fulfilled for the ultrastatic wormhole perturbation. 

At this moment, we are working to include sources of gravitational perturbations. Finally, we can conclude that 
an analysis on the gravitational perturbations to second order, either analytical or numerical, can be more conclusive 
and is on the way to study. 
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